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Abstract 

Based on QCD conformal partial wave expansion to leading order conformal spin accuracy, 
we present the light-cone distribution amplitudes (DAs) of S and A baryons up to twist 
6. It is concluded that fourteen independent DAs are needed to describe the valence 
three-quark states of the baryons at small transverse separations. The nonperturbative 
parameters relevant to the DAs are determined within the framework of QCD sum rule 
method. With the obtained DAs, a simple investigation on the electromagnetic form 
factors of these baryons are given. The magnetic moments of the baryons are estimated 
by fitting the magnetic form factor with the dipole formula. 
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1 Introduction 



Theory of hard exclusive processes in QCD has been studied extensively for several decades. 
The investigation for these transitions provides unmatched opportunities to understand 
the hadron structure, and the theoretical method for the calculation of these processes was 
developed early in 1970's [HE]- In the model of the hard exclusive process, the concept 
of distribution amplitudes (DAs), which are the fundamental nonperturbative functions 
describing the hadronic structure, was introduced. The DAs, physically speaking, describe 
the decomposition of the hadron momentum in parton configurations, which is important 
to make the QCD description of hard exclusive reactions quantitative. 

DAs of mesons have been investigated extensively in the past, some of which have 
been done to high twist accuracy [31 HJ El El [7] . However, the corresponding studies on 
baryons received less attention due to their relatively complex structure, and the existing 
investigations were mainly focused on the nucleon (See [8] for a review). DAs of the nucleon 
and other octet baryons were firstly calculated within the QCD sum rule framework on 
the moments in Ref. [9j. A systematic study on the nucleon DAs were provided in Ref. 
[To] , in which the DAs of the nucleon are investigated up to twist 6. In the paper [TT] . 
the axial vector higher twist DAs of the A baryon were given to leading order conformal 
spin accuracy. Recently Ref. [12] gave DAs of A^ in the heavy quark limit, and Ref. 
[T3] offered a complete analysis of the one-loop renormalization of twist-4 operators of 
the nucleon. In the literature [H], the author presented a description on DAs of helicity 
A = 3/2 baryons with a new approach. Actually, the investigation for DAs of S and A 
baryons was firstly done by Chernyak et al. in Ref. [9j, in which they calculated the DAs 
up to leading twist order only. Nevertheless, more detailed description of the internal 
structures of these baryons needs information on higher twist DAs. The present work is 
devoted to give an investigation on S and A DAs up to twist 6, and then to study their 
electromagnetic form factors as an application. 

Higher twist contributions to DAs come from several physical origins. The first con- 
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tribution is from "bad" components in the wave function and in particular of components 
with "wrong" spin projection. The second one comes from transverse motion of quarks 
in the leading twist components. Finally, higher Fock states with additional gluons or 
quark/ antiquark pairs also contribute to the DAs. It has been known that for mesons, 
contributions due to "bad" components in the quark- antiquark DAs can be described in 
terms of higher Fock states by equations of motion fS^, 'S]. Since quark- ant iquark-gluon 
matrix elements between the vacuum and the meson state are numerically small, contribu- 
tions of "bad" components to mesonic DAs are small enough. However, things are different 
for baryons because equations of motion are not sufficient to eliminate the higher-twist 
three-quark system with additional gluons. At the same time, matrix elements of higher 
twist three-quark operators are large compared with the leading one. Thus the first contri- 
bution is assumed to dominate the DAs rather than the other two origins. For this reason 
we only consider contributions coming from "bad" components in the decomposition of 
the Lorentz structure in this paper. 

The usual description of DAs is based on the conformal symmetry of the massless 
QCD Lagrangian for dynamics dominated on the light-cone. DAs with definite twist can 
be expanded by partial wave functions with the specific conformal spin j. The conformal 
spin of a quark is defined as j = {1 + s)/2, where / is the canonical dimension of the quark 
and s is its spin. For a composite particle, contributions of the higher order conformal spin 
j+n {n = 0, 1, 2,...) are given by the leading contribution multiplied by polynomials which 
are orthogonal over the leading weight function. In this paper, DAs of S and A baryons are 
investigated on the conformal partial wave expansion approach. At first glance, the mass 
terms of the s quark break the conformal symmetry of the QCD Lagrangian explicitly, 
and the SU{3) breaking corrections seem difficult to be included. However, this is not a 
problem as argued by Ball et al. in the case of mesons [5]. The transverse wave functions in 
the conformal expansion are dependent on the scale relevant to the process. If the s quark 
mass is smaller than the QCD scale, the transverse-momentum dependence is not affected 
by the quark mass. Therefore the conformal expansion of the DAs can be carried out 
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safely. In this work, we utilize the method proposed in Ref. [TO] and consider the SU{3) 
flavor symmetry breaking corrections. In fact, the SU{3) flavor symmetry breaking results 
in two different effects, which are isospin symmetry breaking and the corrections to the 
nonperturbative parameters. Unlike the case of nucleon, in which the isospin symmetry 
leads to symmetric relationships to reduce the number of the independent DAs to 8, in our 
case the description of A and S baryons needs fourteen independent DAs that are expanded 
in operators with increasing conformal spin. With equations of motion, the parameters 
of the conformal expansion are expressed in terms of local nonperturbative parameters, 
which need to be determined by nonperturbative QCD methods. In the calculation, we 
expand the DAs to leading order conformal spin accuracy, and use QCD sum rules to 
determine the nonperturbative parameters. 

DAs provide large opportunities to investigate processes connected with the baryons in 
the framework of light-cone sum rule (LCSR) since they are fundamental input parameters 
in this framework [HI [151 [13 [13 [HI [IS]- LCSR is a developed nonperturbative QCD 
method that includes the traditional a la SVZ sum rule [20] technique and the theory of 
hard exclusive processes. The main idea of LCSR is to expand the products of currents 
near the light-cone, and the nonperturbative effects are described by DAs rather than 
condensates in the traditional QCD sum rule [211 [221 [23]. As a simple application, the 
electromagnetic (EM) form factors of the baryons are examined in LCSR with the obtained 
DAs. EM form factors are fundamental objects for understanding the inner structure of 
the hadron. As there are no experimental data available on S and A EM form factors, 
it is instructive and necessary to give an investigation theoretically. In the experimental 
point of view, the EM form factors can be described by the electric and magnetic Sachs 
form factors, and the magnetic Sachs form factor at zero momentum transfer defines the 
magnetic moment of the baryon. It is assumed that the dependence of the magnetic form 
factor on the momentum transfer can be expressed by the dipole formula, therefore after 
fitting the magnetic form factor by the dipole formula, we give estimations on the magnetic 
moments of the S baryons. 
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The paper is organized as follows. Section [2] presents notations and definitions of the 
baryon DAs up to twist 6 by the matrix element of the three-quark operator between the 
hadron state and the vacuum. This section also gives properties of DAs from the symmetry 
relationships. Section [3] gives the leading-order conformal expansion of DAs based on the 
conformal invariance of the Lagrangian. The DAs are simplified to the nonperturbative 
parameters which can be calculated in the framework of the QCD sum rule. Section H] is 
devoted to derive the QCD sum rules for the nonperturbative parameters related to the 
DAs and then present the numerical analysis. Section [5] is a simple application of the 
obtained DAs to investigate the EM form factors of the baryons. We give in this section 
the dependence of the EM form factors on the momentum transfer. After fitting the results 
by the dipole formula, the magnetic moments of the baryons are estimated numerically. 
Summary and conclusion are given in section [61 

2 Definitions of the light-cone distribution amphtudes 
2.1 General classification 

In the SU{3) flavor symmetry limit, light-cone distribution amplitudes of octet baryons 
with quantum number = can be expressed in terms of the matrix element of the 
gauge-invariant operators sandwiched between the vacuum and the baryon state: 

{0\e^^%l{a,z)q2^{a2z)qs';{asz)\X{P)), (1) 

where letters a, P, 7 refer to Lorentz indices and i, j, k refer to color indices, denote the 
quark fields, 2; is a light-like vector which satisfies z"^ = and are real numbers denoting 
coordinates of valence quarks. 

In view of the Lorentz covariance, spin and parity of the baryons, the general decom- 
position of the matrix element is written as [ini [H] : 

m^''%aMQ2jsMq4ia,z)\XiP)) 
= S,MC^p (75X)^ + S2M^C^^ (/75X)^ + V^M {^.C)^^, X, + V2M'' {^.C)^^ (/X)^ 
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+V4M2 (75X)^ + V5M2 (7^0^^ (2a^^^.75X)^ + VgM^ (/75X)^ 

+^1 (f 75^^),^ X, + A2M {ri,C)^p (A), + A,M (7m75C),^ (7^X)^ 

+ (/75C),^ X, + (7m75C),^ (^a^'^z.X)^ + A,M^ (/TsC)^^ (/X)^ 

+%M {a,,C)^^ {a^'l.X)^ + %M (PV^,C)„^ {^"'z.^.X)^ 

+%M' {z-ia^,C)^p (7^75^)^ + {z^ P'' la ,,C) {^^,X)^ 

+%M^ {<y,,C)^^ {a^^ /75X)^ + TgM^ (o^'^^^.75X)^ , (2) 

where X^ is the spinor of the baryon, C is the charge conjugation matrix and a^y = 
'^[ifj.^lu]- In the present investigations denotes baryons with the quantum number 
I{J^) = for T,^ and I{J^) = 0(|^) for A (/ is the isospin, J is the total angular 

momentum and P is the parity). All the functions Si, Vi, Ai, Vj and % depend on the 
scalar product P ■ z. 

The "calligraphic" invariant functions in Eq. ([2]) do not have a definite twist however, 
thus the twist classification needs to be carried out in the infinite momentum frame. Here 
we introduce the second auxiliary light-like vector: 

Vi, = Pi,- J^M— - ' = , (3) 

1 -p ■ z 

so that P — i> p if the baryon mass can be neglected. In the infinite momentum frame, the 
baryon is assumed to move in the positive direction, hence and z~ are the only non- 
vanishing components of p and z. In this frame, terms on twist can be classified in powers 
of p"*", and the baryon spinor is decomposed into "large" and "small" components X'!^ 



and X^ : 
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where the projection operators 

project the spinor onto the "plus" and "minus" components. From the Dirac equation 
f^X(P^ = MX{P), we get the following useful relations: 

^X(P) = MX+(P), A(P) = ^X-(P). (6) 



Using the explicit expressions for X{P), it is easy to see that A~^X = X^ ~ y^jr^ while 
A~X = X~ ~ l/^JjF. By the twist counts in terms of the definition of light-cone 
DAs with a definite twist is given as 

4(0|e^^Vl(ai^)g2^^(«2^)g3^(«3^)|X(P)) 
= S^MC^p (75^+)^ + S^MC^p b^X-)^ + PiM (75C),^X+ + P^M {l,C)^pX- 

+ ¥C)aP (75X+)^ + V2 map (75X-)^ + {^^C)^^ (7^5^+)^ 

+ yM (7xC^)„, (7^5X")^ + 1^5^ (/C^)., (75X+)^ + —n (75X-)^ 

+ f M (7x75C)„, (7^X-)^ + A5— (/75C)„,X+ + —A, {^^,C)^,X~ 

+T, {la^.C)^^ (7^5^+)^ + T, {la^.C)^^ (7^5^-)^ + Ta^ (75X+)^ 

+r4- (75X-)^ + T5— (.ax.C)^, (7^5X+)^ 

+M^K^,C)„,(a^^'75X-)^, (7) 

where an obvious notation = a^^p^Zy, etc., is used as a shorthand and ± stands for the 
projection transverse to z,p, e.g. 7_l7"^ = l^gjluY with 51^^ = ^(^j, - {p^^z^ + z^py)/pz. The 
DAs P = 5*4, Pj, l^, Aj, Tj with a definite twist are classified in Table [TJ Each distribution 



amplitude Fi can be represented as 

F{a^p ■z) = j Vxe-'P'^^''^'''F{xi) , (8) 

where the dimensionless variables Xi, which satisfy the relation < Xj < 1, Yli^i ~ 
1, correspond to the longitudinal momentum fractions carried by the quarks inside the 
baryon. The integration measure is defined as 



J Vx = J dxidx2dxz5{xi + X2 + X3 — 1) . 



(9) 



Comparing Eq. ([2]) and ([7]), the invariant functions Si,Vi,Vi, Ai,% can be expressed in 
terms of the DAs Si, Pi.Vi, Ai, Ti with a definite twist. 

For scalar and pseudo-scalar distributions the following relations hold: 



Si = Si, 2p ■ z S2 = Si ~ S2 , 

Vi=Pu 2p-zV2 = P2-Pu 



(10) 



for vector distributions: 



Vi = Vi, 2p-zV2 = Vi-V2-V3, 

2V3 = V3, 4p- zVi = -2Vi + V3 + V^ + 2V5, (H) 

4p ■ zV, = V^-Vs, {2p ■ zfVe = -Vi + V2 + V3 + V, + V, - , 

for axial vector distributions: 

Ai = Ai, 2p- ZA2 = -Ai + A2-A3, 

2A3 = A3 , 4p ■ zAi = -2Ai - A3 - + 2A5 , (12) 

Ap ■ ZA5 = A3 - A4 , {2p ■ zfAe = Ai - A2 + A3 + A4 - A5 + Ag , 

and, finally, for tensor distributions: 

Ti = Ti , 2p- ZT2 = Ti+T2- 2T3 , 

2T3 = Tj, 2p-z% = Ti-T2-2T7, 

(13) 

2p ■ z% = -Ti + T5 + 2T8 , {2p ■ zf% = 2T2 - 2T3 - 2T4 + 2T5 + 2r7 + 2Ts , 

Ap ■ z% = Tr-Ts, {2p ■ zY% = -Ti + T2 + T5 - Tg + 2T7 + 2Ts . 
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For the identity of the two u{d) quarks in the baryon gives symmetry properties 

of the DAs. The Lorentz decomposition on the 7-matrix structure imphes that the vector 
and tensor DAs are symmetric, whereas the scalar, pseudo-scalar and axial-vector DAs 
are antisymmetric under the interchange of the two u{d) quarks: 

1/,(1,2,3)= 1^^(2,1,3), Ti(l,2,3)= T,(2,l,3), 

5,(1, 2, 3) = -5,(2, 1, 3) , P,(l, 2, 3) = -P(2, 1, 3) , 

A,(l,2,3) = -A(2,l,3). (14) 

The "calligraphic" structures in Eq. ([2]) have the similar relationships. 

For A, the isospin symmetry leads to similar relationships that the vector and tensor 
DAs are antisymmetric, while the scalar, pseudoscalar and axial-vector DAs are symmetric: 

F,(l, 2, 3) = -K,(2, 1,3), T,(l, 2, 3) = -T,(2, 1, 3) , 

5,(1,2,3)= 5,(2,1,3), P,(l,2,3)= P(2,l,3), 

A,(l,2,3)= A(2,l,3). (15) 

2.2 Representation in terms of chiral fields 

This subsection gives the DAs' representation in terms of chiral fields. The discussion is 
mainly about S+ baryon, and the counterparts of the others are similar. In terms of quark 
fields with definite chirality: 

gT(i) = ^(l±75)g, (16) 

the DAs can be interpreted transparently. Projection on the state where the spins of the 
two u quarks are antiparallel, that is u^u^, singles out vector and axial vector amplitudes, 
while the two u quarks are parallel, u^u^ and u^u^, singles out scalar, pseudo-scalar and 
tensor structures. Similar as expressions in Ref. [10], the DAs can be defined in terms of 
chiral fields. The leading twist-3 distribution amplitude can be defined as: 

{0\e'^'(u]{aiz)CMj{a2z)) /sl{asz)\P) = -]^pz ^T.^^ j Vx e-'P'^'^''^ $3(2;,), 
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(17) 

and the distributions for twist-4 are: 

{{)\e'^''(u\{aiz)C ^^u]{a2z)) ^sl{a:,z)\P) = ^pzM jvx e-'^'^'^^^ E^{x,), 

(18) 

and the distributions for twist-5 are similarly written as: 

{0\e''''(uliaiz)Cfh±Mjia2z)'^^^fisiia3z)\P) = -pzM j Dxe-'^^^'^'^' '^^{x^), 

{<d\e'^''(u\{a,z)C^^u]{a2z)^ ^sl{a^z)\P) = ^pzM jvx e-'P'^"^^^^ E^{x,), 

(19) 

and finally the twist-6 one can be expressed as: 

{Ole'^" (ul{aiz)C^u\{a2z)') ^sl{a3z)\P) = --^AP jvx e-'^'^^^'^^ $6(x,). 

(20) 

Due to SU (3) flavor symmetry breaking effects, there are not similar relationships as 
that in Ref. ^I)] from the isospin symmetry to reduce the number of the independent DAs. 
So the following additional chiral fields representations are needed to get all the DAs: 

{<d\e'^^{u\{a^z)C ^MM^^)) HM\P) = \pzM Jvx e-'P'^''^'''{Si - Pi + T3 + Ty) 
{0\e''''(uliaiz)C^flu\ia2z)^ ^slia-iz)\P) = ^pzM ^J:-^^ jvx e-'P'^'="''{S2 - P2 - n + n 



(21) 



and 

{0\e''''(u]{aiz)Cta^,u]{a2z)') 7^/4(MI^) = -2pz ^^^^ Jvx e-'^'^'^^''' Ti{x,), 
{0\e'^^(ulia,z)Cta^,u\ia2z)'^ 7^ ^4(a3^)|P) = -2pz ^^^^ jvx e-'^'^''^''^ T^i^Xi), 
{Q\e''^{u\{a^z)Cia^y^{a2z)) 7^ HM\P) = -M^ ^^^^ jvx 6"'^'^^^^^ T,{x,), 
{0\e'^^(ulia^z)Cta^pu}ia2z)'^ 7^ ^4(a3^)|P) = -M^ Jvx e-'P'^''^''^ nixi). 

(22) 

The twist classification of these additional DAs are shown in Table [H The following 
denotations (^i — Pi + T3 + Ti){xi) = E'^{xi), and (5*2 — P2 — T4 + Ts){xi) = '^'^{xi) are 
adopted for convenience. 

The similar relationships hold for the A baryon under the exchange u, u, s to u, d, s, 
and for the baryon under the exchange u, u, s to d, d, s. 

3 Conformal expansion 

The spirit of the conformal expansion of distribution amplitudes is similar to the partial 
wave expansion of a wave function in quantum mechanics. The idea is to use the conformal 
symmetry of the massless QCD Lagrangian to study the DAs, which allows to separate 
longitudinal degrees of freedom from transverse ones [3l [71 [10] . The transverse coordinates 
are replaced by the renormalization scale, which is determined by the renormalization 
group. The dependence on longitudinal momentum fractions, which is living on the light 
cone, is taken into account by a set of orthogonal polynomials that form an irreducible 
representation of the collinear subgroup SL{2, R) of the conformal group. As to leading 
logarithmic accuracy, the renormalization group equations are driven by tree-level counter 
terms, they have the conformal symmetry. This leads to the fact that components of the 
DAs with different conformal spin do not mix under renormalization to this accuracy. 

The SL{2,R) group is governed by four generators P+, M_+, D and K_, where the 
definitions are used for a vector A: = Af^z'^ and A^ = A^p^ /p - z. The four generators 
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P^, K^, D and M^^, are the translation, special conformal transformation, dilation and 
Lorentz generators, respectively. The generators of the collinear subgroup SL{2, R) can 
be described by the following four operators: 

L+ = -iP+ , L_ = , Lo = -^(D - M_+) , E = «(D + M_+) . (23) 

For a field living on the light cone ^{z), the acting of the above generators on it yields 
the following relations: 

[L^, $(^)] = j(j - miz) , [E, <|.(z)] = (/ - sMz) , 

[E,L2] = 0, [E,Lo] = 0, (24) 

and 

L2 = - Lo + L+L_ , (25) 

where j = {I + s)/2 is called the conformal spin and t = I — s is the twist. In above 
notations, I is the canonical dimension of the quark field and s is the quark spin projection 
on the light-cone. The role of the generator E is analogous to the Hamiltonian in quantum 
mechanics, and the twist corresponds to the eigenvalue of the Hamiltonian. For a given 
twist distribution amplitude, it can be expanded by the conformal partial wave functions 
that are the eigenstates of and Lq. 

For multi-quark states, we need to deal with the problem of summation of conformal 
spins, and here the group is non-compact. The distribution amplitude with the lowest 
conformal spin jmin = ji + 32 + js of a three-quark state is [31 H] 

y>a4xi,x2,x3j r[2^-jr[2j2]r[2j3] ' ' ' ■ ^ ^ 

Contributions with higher conformal spin j = jmm + n {n = 1,2,...) are given by 
multiplied by polynomials that are orthogonal over the weight function fl26|) . In this 
paper, we just consider DAs to leading order conformal spin accuracy. For DAs in Table 
[H we give their conformal expansion: 

<l>s{xi) = 120xiX2Xs4{fi), T^{xi) = 1200:1X2X30:5° (/i) , (27) 
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for twist 3 and 



<l>4(Xi) = 24xiX204(/i) 



T2(Xi) = 24XiX204(/i) 



E'^ixi) = 24x2X3^;° (/i) 



(2^ 



for twist 4 and 



$5(Xi) = 6X30^ (/i) 
S5(Xi) = 6Xi^5(/i) , 
T5(Xi) = 6X305 (/i) 



(29) 



for twist 5 and 



$6(xO = 20°(/x) 



T6(xO = 20;;(/i) 



(30) 



for twist 6. There are altogether 14 parameters which can be determined by the equations 
of motion. 

3.1 DAs of the S baryon 

The normahzation of the DAs of S"*" are determined by matrix element of the local three- 
quark operator. The Lorentz decomposition of the matrix element can be expressed ex- 
plicitly as follows: 

me'''<iOH{0)s'^m^^P)) = V°(/'C)„,(75S+), + V3°(7mC).M7.75S+). 



+Ti°(P'^^a^,C),^(7^5S+)^ + T3°M(a^,C),^(a^'^75S^ 



(31) 



Similar to definitions in Ref. |[10j, the above four parameters can be expressed by the 
following matrix elements: 

(Ole^^'^ [u\0)C A^(0)] 75 ^s'mP) = fMp ■ z) ^E+iP) , 
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(Ole^^'^ [m^(0)C7X(0)] 757^s^(0)|P) = AiMS+(P) , 

(Ole^^'^ [u\Q)Ccj,,u^{id)\ 75fx'^'^s^(0)|P) = A2MS+(P) , 

(0|e*^'^ [MXO)«gV^,M^(0)] 757;.s^(0)|P) = A3M #]+(P) . (32) 

As mentioned above, there are no relations derived from isospin symmetry, so four but 
not three matrix elements are needed to determine the four parameters V^, V3 , T^'', and 
7^°. After a simple calculation, we arrive at the following expressions of , V3 , 7^°, and 
with the four parameters defined in Eqs. fl32l) : 

V? = /s+, V3° = ^(/s+-Ai), 

7'i° = ^(4A3-A2), r3° = l(2A3-A2). (33) 

At the same time, coefficients of operators in Eqs. (l27|) - (!30|l can be expressed to leading 
order conformal spin accuracy as 

<Pl = <Pl = + Ai), = = -^5° = ^(4A3 - A2) , 

= a = J(8A3 - 3A2), <pf = -a' = , 

D 

e;° = ^(12A3-5A2). (34) 
3.2 DAs of the A baryon 

The Lorentz decomposition of the local matrix element of A can be expressed explicitly 
as follows: 

4(0|e^^VjO)4(0)4(0)|A(P)) = 5°a/3(75A)7 + V',{^,C)^^A, 

+^?(/'75C)„/3A^ + ^°M(7^75C)„/3(7^A)^ . (35) 
In order to get the above parameters, the following matrix elements are used: 

(0|e*^''= [u\0)C^, ^d^iO)] ^s'mP) = U{p ■ z) /A(P) , 
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{Q\e'^^ [M'(0)C757^f^'(0)] i^s\Q)\P) = \iMA{P) , 
{0\e'^' [m'(0)C75C^^(0)] s'{0)\P) = X^MAiP) , 

(0|e*^'= [u\0)Cd^iO)] -f5s\0)\P) = \^M^A{P) . (36) 
A simple calculation leads to the following relationships: 

Vl = A2, si = A3. (37) 

To leading order of the conformal spin expansion, coefficients of operators in Eqs. 
(|27|) - (|5n]) for A can be expressed as 



6^ = 0^ = -/a, 0^ = 0^ = -^(/a + Ai), 



^4 = ^5 = \Uk - Al), ^4° = ^5° = A2 + A3, 



a = = A3 - A2. (38) 

4 Determination of the parameters in QCD sum rules 
4.1 QCD sum rules for S baryon 

Determination of the nonperturbative parameters /s, Ai, A2 and A3 can be done in the 
framework of QCD sum rule. The method is carried out from the following correlation 
functions for E"*": 

n.(g') = j ciW^-^(0|T{j,(x)i,(0)}|0), (39) 
with the definitions of the currents: 

ji(x) = e'^^[u\x)C /^^■(x)]75 ^s\x), (40) 
j^{x) = e'^^[u\x)C^,u^{x)]^,^^s\x), (41) 
j3(x) = e'^^[u\x)Ca,,u^{x)]^,a>''s\x), (42) 



14 



and the forth current: 

Ux) = e'^'[u\x)Ciq''a,,u^{x)]^,rs\x). (43) 

Inserting the complete set of states with the same quantum numbers as those of E+, the 
hadronic representations of the correlation functions are given as 

- J so s-q^ 



M - r J so s-r 

2, 2,,2,2 4+M , r P^is) 



U4q') = q'M'Xi^-^^ + / ^ds. (44) 

^ -r J so ^-T 

On the operator product expansion (OPE) side, condensates up to dimension 6 are taken 
into account. To give the sum rules, we utilize the dispersion relationship and assume the 
quark-hadron duality. After taking Borel transformation on both sides of the hadronic 
representation and QCD expansion, and equating the two sides, the final sum rules are 
given as follows: 

4(27r)V|+e"^ = - H s{l - x)'e'^ds - - H ~ ^^^^ ~ e'^dg, (45) 



5 J ml ^ J m 



2 



and 



4(27r)^A?M2e ^ = / ^^[(l - x)(l + x)(l - 8x + x^) - 12x^ lnx]e ^ ds 



b so -773- , 8 ml mlml 

+ - / il-xfe ^ds + -a^(l \ ^ 

&Jmi ' 3 ^ 2M| 2M| 

™4™4 m? pso 



16M| 



_ rso _ s 

)e ^ -2asms / e ^ds, (46) 



and 



2{2nyXlM^e-^'/^B ^ _ / g^-l + 8x- 8x^ + x^ + 12x^ ]nx)e ^ds 

J ml 

+- {1 - x){4: - 7x)e ^ ds + Urustts e^ds, (47) 

3 J mi J mi 
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and 



/•so 

(47r)^A2MV^^'/^^s = / sm{l-x){l+x){l-8x + x'')-l2xHnx] 

J ml 

+-(l-x)^}e AT^fisH / (1 -x)(ll -5x-4x^)e '^cis 

5 12 



2 2 2 A A 

2, rriQ mQrrig mQ-ml 



2 



+ 16a^( 1 H )e 

/•so _ s 

-8m,a, / e "J(is. (48) 

where the notation x = m1/s is used, and the other parameters employed are the stan- 
dard values: a = -{27ry{uu) = 0.55 GeV^ b = {2Tif{a,G'^/^) = OAT GeV^ a, = 
— (27r)^(ss) = 0.8a, and {ugc<y ■ Gu) = 0.8{uu). 

As the usual way of the sum rules, the auxiliary Borel parameter should have a 
proper range in which the results of the sum rules vary mildly with it. On the one hand the 
Borel parameter is expected to be large so that the higher order dimensional contributions 
are suppressed, and on the other hand the Borel parameter needs to be small enough to 
suppress the higher resonance contributions. Fig. [1] shows the dependence of the above 
parameters on the Borel parameter M^. The window of Borel parameter is choose as 
1 GeV^ < < 2 GeV^, in which our results are acceptable. 

To determine the relative sign of /s+ and Ai, we give the sum rule of /s+A^: 

- ^■'^ m /"*° - " 

4(27r)Vs+AtMe ^ = — / s[{l - x){3 + 13x - 5x^ + x^) + 12x\nx]e ^ds 

3 Jml 

Am, r 1(1 -x)[l+ (l-^)(2-5x) j^--^^^ _ 4^^ ro e'^ds. (49) 
6 J^2 s 3a; 3 J^2 

It is similar that the relative sign of A2 and A3 to Ai can be given by the following two 

sum rules: 

(27r)^(AiA; + AiA2)M^e ^ = -12m,a / {I + ml/ s){l - xfe ^ds, (50) 



and 



fSo 

{2tt)\XiX; + AtA3)MV^ = - / {as{l - xy{2 + x) 



m 



2 
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TTl Oj 3 — ^ 

+^[1 - -(1 - x)(l + x) + (1 - x)(13 - 25x + 2x^)]}e ^ds. (51) 

Fig. [2] gives the dependence of the above sum rules on the Borel parameter M^. 
The final numerical values of the coupling constants of S are: 

/s = (9.4 ± 0.4) X 10-2 GeV^ Ai = -(2.5 ± 0.1) x 10"^ GeV^ 

Aa = (4.4 ±0.1) X 10^2 GeV^ A3 = (2.0 ± 0.1) x 10^^ GeV^ (52) 

In Ref. [9], Chernyak et al. have calculated nonperturbative parameters relevant to 
leading twist contributions. Their numerical results are — 0.51 x 10~^ GeV^ and 
l/^l ~ 0.49 X 10^^ GeV^. By contrast, we use Eqs. (!32|) and (1331) to give our estimations: 
|/s| = 0.94 X 10-2 GeV^ and |/|| = 0.60 x 10"^ GeVl In comparison with their results, 
our predictions are lager in absolute values. 

4.2 QCD sum rules for A baryon 

The sum rules of the A baryon parameters begin with the following correlation functions: 

liM') = I d'xe''-^0\T{Ux)m}\0), (53) 
with the definitions of the currents: 



Jll 


[x) 


= e'^'' [u'\ 


[x)C^,^d\x)] ^s'ix), 


(54) 


J2I 


[x) 


= e'^^ [u'\ 


[x)C^,7,d^ix)] Ys'ix), 


(55) 


J3I 




= e'^^ [u'\ 


[x)Cj,d^{x)] s\x), 


(56) 


J4I 


[x) 


= e'^'' [u'\ 


[x)Cd^{x)] ^,s'{x). 


(57) 



The similar processes as in the above subsection lead to the following results: 

(47r)Vle ^ = - / s{l-xfe^ds-- -x{l - x){l - 2x)e ^ ds, (58) 

5 J ml 3 s 

and 

A{2TT)^XlM^e ^ = - / s^[{l-x){l + x){l-8x + x^)-12x'^\nx]e ^ ds 

2 J ml 
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h , -773- , 4 ml mlml 

+— / il-xfe ^ds--a^(l ° ° ' 



12 ^ 3 ^ 2M| 2M^ 

■-^)e '^-m.a, / e "Wrfs, (59) 



and 



_ /-SO _ s 

A{ATT)^XlM^e ^ = / s2[(l-x)(l+x)(l-8x + x2) - 12x2lnx]e "Wds 



+ - / (1 -x)(l + 5a;)e ^ ds + —a\l - ° ° ' 



3 7^/ ' 3 ^ 2M| 2M| 

-^^)e "W-4m,a, / e ^rfs, (60) 



and 



_ rsQ _ s 

4(47r)^A^M2e "I = / s^[(l - x)(l + x)(l - 8x + x^) - 12^^ lnx]e "^c/s 
+ - / (l-x)(l + 5x)e ^ds + —a^{l 

3 im? 3 



b _ , --rfr , 32 2^-, "^0 mlml 

2Ml ~ 2M% 



+ TKU^^^ ^-4m,a, / e ^rfs, (61) 
16M« J^2 

and the sum rule of /a A]; is 

(47r)VAAlMe ^ = -m, / s[(l - x)(3 + 13x - Sx^ + x^) + 12xlnx]e '^ds 

b pi. xr (1 -a;)(2-5a;), --fr , 8 p --f^ , , , 
+ -m, / -(l-x)[l + ^ '-]e^ds+-as e^ds. (62) 

Note that the sum rules of A2 and A3 are the same. In Fig. [3l the sum rules of the 
parameters on the Borel parameter are shown. The final numerical results for the 
parameters of A are: 

/a = (6.0 ± 0.3) X 10-^ GeV^ Ai = (1.0 ± 0.3) x 10"^ GeV^ 

IA2I = (0.83 ± 0.05) X 10~2 GeV^, IA3I = (0.83 ± 0.05) x lO'^ GeV^. (63) 

In the above results, /a and Ai have the same sign, which is different from that shown in 
Ref. [H]. The relative sign of Ai and A2 cannot be determined by the method presented 
above. Here we only list the absolute values of the two parameters. 
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The numerical estimations by Chernyak et al. are: |/a| — 0.63 x 10 ^ GeV^ and ^ 
0.063 X 10-2 GeVl Our result on /a is /a = 0.60 x 10"^ GeV^ but in our calculations /J is 
zero, which is different from that of Chernyak. The deviation is due to the tensor structure 
of the baryon, which disappears at the leading order of the conformal expansion because of 
the isospin symmetry. In our approach, the determination of the tensor coupling constant 
relies on the next conformal expansion of the DAs, which corresponds to the sum rules of 
higher order moments in Ref. [9]. 

5 Applicationrelectromagnetic form factors of the baryons 
with hght-cone QCD sum rules 

5.1 LCSR for the electromagnetic form factors 

The EM form factors of hadrons are fundamental objects for understanding their internal 
structures. There were a lot of investigations on various hadrons both experimentally and 
theoretically, including meson [Ml EHl ESI EH EHl EH] and baryon [30l EH [32l ESI El ESI 
[36l EH EH El] . While as there were few experimental data and theoretical investigations, 
the EM form factors of the baryons, such as S and A and so on, have not received much 
attention in the past years. Chiral perturbation theory and the chiral quark/soliton model 
have been used to study the baryon EM form factors at low momentum transfer ^01 UT] . 
T. Van Cauteren et al. have investigated the electric and magnetic form factors of these 
baryons in the relativistic constituent quark model [42j. In the previous work [43j we gave 
an investigation on the EM form factors of the A baryon. In this section, the EM form 
factors of the E baryon are investigated at moderately large momentum transfer within 
the framework of light-cone QCD sum rule method, and the magnetic moments of the 
same baryons are estimated by comparing our results with the existing dipole formula. 

The matrix element of the electromagnetic current between the baryon states can be 
expressed as the Dirac and Pauli form factors FiiQ"^) and F2{Q'^), respectively: 

{^iP,s)\j;-iOmP',s')) = ^P,s)[^,F^iQ')-^^F,iQ')nP',s'), (64) 
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where j^™' = e„'U7^M+esS7^s is the electromagnetic current relevant to the hadron, and P, s 
and P', s' are the four-momenta and the spins of the initial and the final S baryon states, 
respectively. From the experimental viewpoint, the EM form factors can be expressed by 
the electric Ge{Q^) and magnetic Gm{Q^) Sachs form factors: 



GM{Q')=F^{Q^) + F2iQ 

4M2 



Ge{Q') = F,{Q') - S^F,{Q% (65) 



and at the point = the magnetic Gm{Q'^) form factor gives the magnetic moment of 
the baryon: 

G'm(O) = /is. (66) 

In order to evaluate the magnetic moment of the baryon from its EM form factors, the 
magnetic form factor Gm^Q"^) is assumed to be described by the dipole formula: 

—GMm = = Gum- (67) 

As there is no information about the parameter from experimental data, the two 
parameters rn^ and /is are estimated simultaneously by fitting the magnetic form factor 
with the dipole formula (!67|) . 

The calculation mainly focuses on baryon and it is similar for the calculation of 
S^. The process of the derivation begins with the correlation function: 

T^(P,g) = z j d^xe^'^-^(0|T{js+(0)jr(x)}|S+(P,.)), (68) 

where the interpolating current js+ is chosen as Eq. fHOj) . The hadronic representation 
of the correlation function is acquired by inserting a complete set of states with the same 
quantum numbers as those of S"*": 

+ {P' ■ zF,{Q') + ^P2(Q^))^]S+(P, .) + ..., (69) 
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where P' = P — q, and the dots stand for the higher resonances and continuum con- 
tributions. The correlation function is contracted with to get rid of contributions 
proportional to which is subdominant on the light cone. On the theoretical side, the 
correlation function fl68|) can be calculated to leading order of as 

z;r^^ = (p.^)2(/s+)J4e„ / rf«2 ^{SoM + T m^iM 

Jo s — p [s — p ) 

1 

-2 (g_p.2)2 -^2(a2)} + 2e, das^—^iCoias) 
:Ciias)-2- -—C2ias)} 



^ Jo a2[s-P'^Y 
1 

[s - P'^) Jq a3(s2 - P'^)^ 



0^2 



(«3) 



+2 , ^' c^ia-s)}], (70) 

{S2 - P'^) J 

where s = (1 - a2)M2 + ^i^Q^ ^nd S2 = (1 - a^)M'^ + ^^^Q^ + ^. Here the following 
notations are used for convenience: 

Bo{a2) = / daiVi{ai,a2,l - ai - 02), 
Jo 

B,{a2) = (2Vi -V2-V,-V,- V,){a2), 
B2{a2) = {-V, + V2 + Vs + n + V,- Ve){a2), 

/•l-«3 

Coias) = / daiVi{ai,l - ai - a3,a3), 
Jo 

Ci(«3) = {2Vi -V2-V3-V,- V^){a3), 
C2(a3) = {-Vi + V2 + V3 + V^ + V^- Ki)(a3), 
/^i(a2) = {Vi-V2-V3){a2), 

Eiias) = {V,-V2-V3){a3), (71) 
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where 



Vi{oi2) — / / daiVi(ai, a2,l — Oil — 0.2), 

Jo Jo 

~ ^OL2 f-a'^ r^-<^2 

Vi{a2) = / da2 / da2 / daiVi{ai,a2,l — ai — a2), 
Jo Jo Jo 



'»a3 pl—a'^ 

^i(ct3) = / da'^ / daiVi{ai,l — ai — a'^ja'^), 



JO 



"as pa'^ /"I— Q3 

Vi{as) = I da'^ / da'^ / daiVi{ai,l — ai — a'^^a'^). (72) 



^0 ^0 



Then equating both sides of the Borel transformed version of hadronic and theoretical 
representations with the assumption of quark-hadron duahty, the final sum rules are given 
as follows: 

2/s+i^i(Q')e ^ = 4e„ / da2e ^ 5o(a2) + tf^5i(«2) - T:r^52(«2) 

+2e, / (iaae C^[a^) + — ^Ci(a3) - _ ^Calag) \ 



+2e.e , , ^, :r-C2{aso) 



(73) 



for Fi(g2) and 

h+F2{Q')e~^ = M'Ueu / da^e rrrl - ^1(^2) + TrT^2(a2)} 



' 020 

fO Q;2q 

fo Q;LM^ d ^ , , 0:20 

-4e«e ^-^i^^x-^2(a2o)- 
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£2_ 1 

Jaso asMs Mb 



-C2(«3o) 2 /I #2 I ^/^2 I 2 1 



aioM2 + Q2 + ^2 ^ct3o aioM2 + Q2 + ^2 

(74) 

for F2(g2). 

5.2 Numerical analysis 

In the numerical analysis, the continuum threshold is chosen as Sq = (2.65 — 2.85) GeV^. 
The masses of the S baryons from Ref. P] are Ms+ = 1.189 GeV and M^- = 1.197 GeV. 
The parameters /e and Ai are used as the central values in Eqs. (1521) . For the auxiliary 
Borel parameter M^, there should be a region where the sum rules are almost independent 
of it. To choose a platform for the Borel parameter, we should suppress both resonance 
contributions and the higher twist contributions simultaneously. Fig. H] shows the depen- 
dence of the magnetic form factors on the Borel parameter at different points of Q^. Our 
results are acceptable in the range 2.0 GeV^ < Ml < 4.0 GeVl 

The estimation on the magnetic moment of the baryon comes from the fitting of the 
magnetic form factor by the dipole formula fp7|) . In the following analysis the Borel 
parameter is chosen to be = 3 GeV^. Fig. Ogives the dependence of the magnetic 
form factor G'm(Q^) on the momentum transfer at different points of the threshold Sq. The 
figure shows that Gm{Q^) decreases with Q^, which is in accordance with the assumption 
in Eq. (1671) . To estimate the magnetic moment numerically, the magnetic form factor 
Gm{Q^) is fitted by the formula /is/(l + Q'^/^q)'^, which is described by the dashed lines 
in Fig. O From the figures the magnetic moment of S+ is given as = (3.13±0.10)/i7v, 
and the estimation of the other parameter is ml = (0.86 ± 0.04) GeV^. 

The similar process is carried out for the numerical analysis of S^. The estimation of 
the S~ magnetic moment is shown in Fig. [6], which are the magnetic form factors of the 
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baryon at different threshold and the fittings by the dipole formula. The numerical values 
from the fittings are /is- = -(1.59 ± 0.02)/i^ and = (0.78 ± 0.03) GeVl 

Table [2] lists magnetic moments of the two baryons from various approaches: data 
from Particle Data Group (PDG) [H]; QCD sum rules |15] (SR(1) for x = -3.3 and 
SR(2) for X = —4.5); QCD string approach (QCDSA) ^46j; chiral perturbation theory 
(xPT) HH; Skyrme model (SKRM) [H]; light cone sum rules Hg (LCSR(l) for x = -3.3 
and LCSR(2) for x = —4.5). The table shows that our results are larger in absolute 
values than the others. This may partly lie in the fact that more detailed information 
on DAs calls for higher order conformal expansion. At the same time, the choice of the 
interpolating currents may affect the results to some extent [I5l [16]. The estimation is 
expected to be better if more information about the DAs are known and higher order QCD 
coupling 0{as) effect are included. 

Finally, the Q^-dependence of the physical value Gm/iJ'T.Gd is given in Fig. [71 In 
the numerical analysis, the input parameters used in the dipole formula f[671) are the 
central values obtained above, which are = 0.86 GeV^ for T.'^ and ttiq = 0.78 GeV^ for 
S^, while the magnetic moments used come from Ref. [44], which are = 2.458/iAr and 
/is- = -1.160/XAr. 

6 Summary 

In this paper we present the DAs of baryons with quantum number I{J^) = 1(|''^) (for 
E^) and I{J^) = O(l^) (for A) up to twist 6. We find that fourteen independent DAs 
are needed to describe the structure of the baryons. The method employed is based on 
the conformal partial wave expansion, and the nonlocal nonperturbative parameters are 
determined in the QCD sum rule framework. Our calculation on the conformal expansion 
of the DAs is to leading order conformal spin accuracy. Compared with the previous work 
[TT] , the calculation on the A baryon gives DAs of all other Lorentz structures besides axial- 
like vector structures. Another new result is that the relative sign of the two parameters 
/a and Ai are positive. 
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With the DAs obtained, the EM form factors of S are investigated in the range 
1 GeV^ < < 7 GeV^. We assume that the magnetic form factor can be described 
by the dipole formula. Fitting the result by the dipole formula, the magnetic moments of 
the baryons are estimated, which are = (3.13±0.10)/i7v, and yUs- = — (1.59±0.02)/iAr. 
Compared with values given by Particle Data Group [44j, our results are larger in abso- 
lute values. This shows that the calculation needs more detailed information on the DAs, 
which may come from higher order conformal spin contributions, and at the same time 
the choice of the interpolating currents may also affect our estimation to some extent. 
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Appendices 

In the appendices we give our results on the DAs of S and A explicitly. As the definition 
in Eq. ([7]), our results are listed in the following subsections. 

A DAs of the S baryon 

Twist-3 distribution amplitudes of S are: 



'I 



Ai{x,) =0, 



Ti(Xi) = 12OXiX2X303° . 



(75) 



Twist-4 distribution amplitudes are: 



Si{Xi) = 6(X2 - Xi)x3(^4 + ^4°) , 



=6(x2-Xi)x3(e4-e4°), 



V2{X,) = 24x1X20: 



'4> 



A2{Xi) = 0, 



1^3 (X,) = 12X3(1 -Xs)^; 



'I 



Asixi) 



12X3(Xi - X2)'04 , 
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T2iXi) = 24xiX20l°, 

r7(x,) = 6x3(1 -a;3)(C-e4)- 
Twist-5 distribution amplitudes are: 

VsiXi) = 6X30^, 

Ux,)^-l{x,+X,){^',' + e,): 

Ts{x,) = l{x^ + X2m','-e,). 
Finally twist-6 distribution amplitudes are: 

V,{xi)^2<f>l 

B DAs of the A baryon 

Twist-3 distribution amplitudes of A are: 

V,{xi)^0, 

Ti{xi) = 0. 
Twist-4 distribution amplitudes are: 

S,{xi) =6x,{l-Xs){e,+^T), 

V2{Xi) = , 

V3{Xi) = 12{xi - X2)x3ljjl , 

T2{xi) = 0, 

T^ixi) = -6(a;i-a;2)a;3({4° + e4°)- 



n{xi)^6xsil-xs){eA + C'), 

(76) 

P2ix^) = l{x,-X2mt-^','), 

Ai{xi) = 3{xi - xa)^/'" , 
A5{xi)^0, 

TsiXi) = 6X305°, 

(77) 

Ae{xi) = 0, 

(78) 



Ai{Xi) = -120X1X2X303, 

(79) 

p,(x,) = 6(i-x3)(e4°-e;°), 

A2{Xi) = -24xiX204 , 
A3(Xi) = -12X3(1 -X3)V'4, 
r3(xO = 6(X2 - Xi)x3(-e4 + ^4°) , 

(80) 



Twist-5 distribution amplitudes are: 

V,{x,) =0, 

Ts{x,) ^ -^ix, - x,){e, - ^'s') . 
Finally twist-6 distribution amplitudes are: 

Ve{xi) = Q, 
n{xi) =0. 
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Fig. 1. Dependence of the four parameters /s, Ai, A2 and A3 of E 
on the Borel parameter M^. The hnes correspond to the threshold 
Sq = 2.65 — 2.85 GeV^ from the bottom up. 

Fig. 2. Sum rules of the relative signs of the parameters on the 
Borel parameter. The threshold is used as Sq — 1.66^ GeV^. 
Fig. 3. Dependence of the four parameters /a, Ai, A2 and /aA^ 
of A on the Borel parameter M|. The lines correspond to the 
threshold Sq — 2.45 — 2.65 GeV^ from the bottom up. 
Fig. 4. Dependence of the form factor Gm{Q'^) of S on the Borel 
parameter at different momentum transfer. The lines correspond 
to the points = 1 ^ 2 , 3 , 5 , 6 GeV^ from the up down (left E+) 
and from the bottom up (right E~), respectively. 
Fig. 5. Fittings of the form factor Gm(Q^) by the dipole formula 
/is+/(l + Q^/'tTT'oY- The dashed lines are the fittings, and figures 
(a) , (b) correspond to the threshold Sq — 2.65 , 2.85 GeV^, respec- 
tively. 

Fig. 6. Fittings of the form factor Gm{Q^) by the dipole formula 
l^-i^-fi^ + Q^/n^o)^- The dashed hnes are the fittings, and figures 
(a) , (b) correspond to the threshold Sq — 2.65 , 2.85 GeV^, respec- 
tively. 

Fig. 7. The Q^-dependence of the form factor Gm/(a*sGd)- The 
hnes correspond to the threshold Sq = 2.65 , 2.75 , 2.85 GeV^ from 
the bottom up. The left corresponds to E+ and the right corre- 
sponds to E~. 
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Table 1: Independent baryon distribution amplitudes that enter the expansion in Eqs. 
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Light-cone projection 
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Table 2: Magnetic moments of the E baryons from various models 



/U(/iiv) 


PDG 


SR(1) 
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-1.16 
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-0.8 


-1.1 


-1.59 
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